clear
clc
close all


data.boot.soleLength = 30e-2;
data.binding.model = 'bulldog';
data.binding.pivotPosition = 2;

Fh = .1;

aBs = linspace(-.1,.15,200);
aTs = linspace(-.1,.1,201);

for i = 1:length(aTs)
    for j = 1:length(aBs)
        res = basicBoot([aTs(i),aBs(j)],.1,data);
        
        rt(i,j) = abs(res(1));
        rb(i,j) = abs(res(2));
    end
end

figure
mesh(aBs,aTs,rt)
ylabel('Toe Angle')
xlabel('bellows angle')
zlabel('toe moment')

figure
mesh(aBs,aTs,rb)
xlabel('bellows angle')
ylabel('toe angle')
zlabel('bellows moment')